Evidence for in-medium modification of (f) meson at normal nuclear density 
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Invariant mass spectra of e^e~ pairs have been measured in 12 GeV p + A reactions to detect 
possible in-medium modification of vector mesons. Copper and carbon targets are used to study the 
nuclear-size dependence of e'^e~ invariant mass distributions. A significant excess on the low-mass 
side of the (f) meson peak is observed in the low I3j{— 13/^^1 — 13'^) region of (f> mesons ( < 1.25 
) with copper targets. However, in the high region {0y > 1.25), spectral shapes of (j) mesons 
are well described by the Breit-Wigner shape when experimental effects are considered. Thus, in 
addition to our earlier publications on p/ui modification, this study has experimentally verified 
vector meson mass modification at normal nuclear density. 

PACS numbers: 14.40.Cs, 21.65.+f, 24.85. +p, 25.40.Ve 



The properties of hadrons, which are composite parti- 
cles of quarks and gluons, have been measured and deter- 
mined in the past in many experimental studies. These 
properties, such as mass and decay width, represent vac- 
uum expectation values. Therefore, in principle, they can 
differ when the vacuum itself changes, i.e., when they are 
in hot and/or dense matter. 

A possible and interesting experimental approach to 
the in- medium properties of hadrons is to measure the di- 
leptons from the decays of hadrons in hot/dense nuclear 
matter, in which the chiral symmetry could be restored 
leading to changes in those properties. The present ex- 
periment, KEK-PS E325, was conducted at the KEK 
12-GeV Proton-Synchrotron, in order to search for in- 
medium mass modifications of vector mesons in the re- 
action 12 GeV p-\-A^n,uj,(j)-\-X^e+e--\-X'. In our 
earlier publications [li Q j we reported the mass modifi- 
cation of p and u mesons in a nuclear medium. In the 
present paper, we report new results of the meson. 

In the case of the p and lo mesons, it is very difficult to 
distinguish between matter effects affecting only one of 
the two mesons, or both, due to the overlapping of two 
peaks with different widths [l|, Q ■ The natural width of 
the (j) meson is narrow (4.26 MeV/c^) without other res- 
onances in the vicinity; therefore, we can examine the 
possible mass modification more clearly. When the prop- 
erties are modified by the medium, the observed spec- 
trum contains two components: in-vacuum decays with 
a normal mass distribution and in-medium decays with 
a modified mass distribution. The latter component can 
cause an excess around the (j) meson peak in the invariant 
mass spectra. 



Theoretically, the possibility of the decrease in the 
mass of light vector mesons in a nuclear medium was 
first pointed out by Brown and Rho 3]. Thereafter, 
many theoretical studies were conducted. Hatsuda and 
Lee calculated the density dependence of the (f> meson 
mass based on QCD sum rules According to their 
calculation, the expected mass decrease for the (j> meson 
at normal nuclear density is 20-40 MeV/c^. As for the 
decay width of the (j) meson, some theoretical calculations 
predict the broadening of the width by a factor between 
five or six [H, \§\ and ten at normal nuclear density. 
When the width of the (p meson broadens by a factor of 
ten, the lifetime of a (/) meson, cr^, in a nucleus is reduced 
from 46 to 5 fm and the probability of in-medium decay 
increases, thereby reflecting the in-mcdium properties. 

Although several experimental reports on the in- 
medium modification of p and tu mesons exist, including 
our reports P, 0, [^, [lO, T7| ^ , experimental informa- 
tion on the (f> meson modification is very limited. The 
CERES/NA45 experiment measured dielectron spectra 
in high-energy heavy-ion collisions the HADES ex- 
periment measured dielectron spectra in 2 GeV-A C-C 
collisions [l^]- However, in both experiments, the mass 
resolution and/or statistics for the mesons were not 
sufficient to draw a definite conclusion regarding the 
6 meson mass modification. The NA60 collaboration 



^ Recently, the GSI S236 collaboration reported an enhancement of 
the pion-nucleus potential in a deeply bound pionic 115>119.123 g^j 
nuclei as a possible signature of chiral symmetry restoration [l2j . 
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measured the dimuon spectra in 158 GeV-A In-In col- 
lisions ll|, and the PHENIX experiment at RHIC re- 
spectrum in Au-|-Au collisions at 
Recently, the LEPS coUabora- 



ported the 

= 200 GeV [l^ 
tion reported a possible ct^at modification in medium by 
measuring the A-dependence of the </> photo-production 
yields in the K^K~ decay mode [15|- Thus far, no clear 
evidence for the modification of the </) meson mass has 
been observed in the above experiments. The result de- 
scribed in the present paper is the first positive signal of 
the meson modification. 

Detector elements relevant to our analysis are briefly 
described as follows. For further details of the E325 spec- 
trometer, see |16j. It comprises two arms with electron ID 
counters and kaon ID counters that share a dipole magnet 
and tracking devices. The typical acceptance in the lab- 
oratory frame was 0.5 < rapidity < 2.0 and 1 < /37 < 3 
for e'^e~ pairs. In the present paper, we report analysis 
results with e+e^ -triggered data collected in 2001 and 
2002. A primary proton beam with a typical intensity of 
9 (7) X 10* per 1.8-sec spill in 2001 (2002) was delivered 
to targets located at the center of the magnet. In order 
to observe the nucleus-size dependence, we accumulated 
data by using two types of targets, carbon and copper. 
In 2001, one carbon and two copper targets were used 
simultaneously, while in 2002, one carbon and four cop- 
per targets were used simultaneously. The thickness of 
each copper target was 73 mg/cm^ and that of the car- 
bon target was 92 (184) mg/cm^ in 2001 (2002). They 
were aligned along the beam axis and separated typically 
by 46 (23) mm in 2001 (2002). 

In order to reproduce the observed invariant mass spec- 
tra, we performed a detailed detector simulation using 
Geant4 



17|. 



All the experimental effects that affect 
the invariant mass spectrum, such as multiple scattering 
and energy loss including the external Bremsstrahlung 
of particles, tracking performance with chamber resolu- 
tion, and misalignment of tracking devices, were consid- 
ered. The effect of internal radiative corrections was also 
taken into account according to [l3| . The mass resolution 
of ^ e+e~ was estimated to be 10.7 MeV/c^. 

We reconstructed the masses of the mesons from the 
measured momenta of the and e~. Figure [T] shows 
the obtained invariant mass distributions. We divided 
the data into three parts based on the (]j values of the 
observed e+e~ pairs, (3j < 1.25, 1.25 < < 1.75, and 
1.75 < We fitted each mass spectrum with a reso- 
nance shape of — > e'^e^ and a quadratic background 
curve. For the (p meson resonance shape, we used the 
Breit-Wigner curve M^{m) oc l/((m - ruo)^ + (ro/2)^) 
with pole mass mo = 1019.456 MeV/c^ and decay width 
Fq — 4.26 MeV/c^ convoluted over the detector response 
in the simulation according to the kinematical distribu- 
tions of the (j) mesons in each f3j region. The kinemat- 
ical distributions of the (f> meson were obtained by the 
nuclear cascade code JAM [l^, which reproduced well 
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FIG. 1: Obtained e'^e" distributions with the fit results. The 
target and /37-region are shown in each panel. The points with 
error bars represent the data. The solid lines represent the 
fit results with an expected (j> e^e~ shape and a quadratic 
background. The dashed lines represent the background. 
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FIG. 2: Kinematical distributions of e^e~ pairs in the mass 
region 0.95 < Mee < 1.05 GeV/c^. Acceptance was not cor- 
rected. Distributions of (a) j3'y and (b) contours in transverse 
momentum and rapidity. In the plot (a), points represent 
data and the line represents simulation result using the nu- 
clear cascade code, JAM [l9l |. 



the observed distributions as shown in Fig. [2l The rela- 
tive abundance of the 4> mesons A^^, and the parameters 
of the quadratic background were obtained from the fit. 
The fit region was from 0.85 to 1.2 GeV/c^. The car- 
bon data were well reproduced by the fit in all the (3^ 
regions. On the other hand, the copper data in the low- 
est (3^ region contradicted the applied resonance shape 
at 99% C.L. due to the visible excess on the low mass 
side of the 4> meson peak. 

To evaluate the amount of excess A^cx, we repeated 



3 



TABLE I: The fit /ndf values (a) including and (b) ex- 
cluding the excess region: 0.95 — 1.01 GeV/c^. 





P'y range 


X^/ndf (a) 


X^/ndf (b) 




< 1.25 


36/50 


31/41 


c 


1.25 - 1.75 


63/50 


54/41 




1.75 < 


46/50 


37/41 




< 1.25 


83/50 


57/41 


Cu 


1.25 - 1.75 


43/50 


38/41 




1.75 < 


55/50 


53/41 



the same fit procedure except that the mass region from 
0.95 to 1.01 GeV/c^, where the excess is observed, was 
excluded from the fit. The obtained x^/ndf' s including 
and excluding the excess region are shown in TABLE [D 
We obtained A'cx by subtracting the integral of the fit re- 
sults from that of the data in the mass region from 0.95 
to 1.01 GeV/c^. To evaluate systematic errors, we varied 
the fit region, binning, mass resolution and mass scale, 
changed the background curve from quadratic to cubic, 
changed the Breit-Wigner shape from the non-relativistic 
one to the relativistic one, and modified the kinematical 
distribution of (j) mesons in the simulation. Then we av- 
eraged the obtained values of A'ox and iV^ to obtain the 
listed values in TABLE |TT1 The systematic errors repre- 
sent the maximum deviations from the averaged values. 
For the statistical errors we selected the largest values in 
above fit conditions. The excess is statistically significant 
for the copper target data in the lowest [3^ bin, whereas 
it is marginal for the carbon target data. This excess is 
considered to be the signal of the mass modification of 
the (j) mesons in a target nucleus because such an effect 
should be visible only for slow (j) mesons produced in a 
large target nucleus. 

In the mass region of the p meson the excess was ob- 
served both in C and Cu targets but in the case of 
the (j) meson we observed the excess only in Cu targets 
data with small (3^. The p-meson lifetime (cr ~ 1.3 fm) 
is shorter than the typical nuclear size and a significant 
portion of p mesons should decay in both C and Cu nu- 
clei. However, the 0-meson lifetime (cr ~ 46 fm) is much 
longer, so even if it is modified in nuclear medium, such 
an effect should only be visible for a larger nucleus and 
for (/)'s with small (3^. 

We attempted to reproduce the observed mass spec- 
tra with a Monte-Carlo-type model calculation that in- 
cludes an in- medium mass modification of the 6 mesons 
based on the theoretical predictions of 0, H, la 0] ■ 
assumed the density dependence of the 4> meson mass 
as m^{p)/mci,{0) = 1 — /ci(p/po), where po is the nor- 
mal nuclear density [^J. To reproduce the large amount 
of excess in our data (22% for slow (f> in the copper 
target data), it is necessary to introduce a broadening 
of the total width of the (j) (rj^°'), or at least of the 
partial width for the (p — > e+e^ decay (F^'^). When 



TABLE II: Numbers of </!> mesons [A'^^] and excess [N^x] eval- 
uated by the fits (see text). Error values shown are statistic 
(first) and systematic (second). 





range 


A^cx A-^ A^ex/(Afcx + A^) 


c 


< 1.25 
1.25 - 1.75 
1.75 < 


6 ±17^6 257±26l? 0.02 ± 0.06l[];oi 
-4±26l"2 547±4lt]_l -0.01 ± 0.05lo:o? 
39±42l25 1076 ±641^5 OM±OMt°°n 


Cu 


< 1.25 
1.25 - 1.75 
1.75 < 


133 ± 28+12 464 ± 38^5 0.22 ± 0.04lo;o^ 
35 ± 33i?2 588 ± 47lg'' 0.06 ± 0.05j;J5;^i 
21±48l 29 1 367 ±72l27 0.02 ± O.OSl^? 



no broadening is introduced, the expected rate of in- 
nucleus decay is just 6% for the <p mesons produced 
in copper nuclei with (3j < 1.25. For the density de- 
pendence of the total width broadening, we assumed 
Tl°'ip)/Tf\0) - 1 + fcr*(p/po); kl°' > 0. For Tf we 
assumed T^{p)/T'^{0) = 1-1- A;|^(p/po), and examined 
following two cases: (i) the branching ratio F^'^/F^"* re- 
mains unchanged in the medium; fcl"^ = kl°*, and (ii) 
doesn't increase in the medium; = 0. It should be 
noted that the ratio A^^l^^^ /^^^^ee increases with F^°* 
even in the case (ii). Here, N^^^ {N^^^J denotes the 
number of in-medium (in-vacuum) (j) — s- e+e~ decays. 
This is because both N^^^ and N^^^^ decrease but the 
latter decreases faster. 

We considered that <j) mesons were generated in the tar- 
get nucleus according to the nuclear density profile. This 
is because we measured the mass-number dependence of 
the (f) meson production cross section as (j{A) cx A"^ [20| . 
Generated cf) mesons were traced until the decay point 
with the modified pole mass and decay width according 
to the nuclear density. The decay probability increases in 
the medium due to the width broadening. We used the 
Woods-Saxon distribution for the nuclear density profile: 
p/po oc (1 + exp((r - R)/t))-^, where R = 4.1(2.3) fm, 
and r = 0.50(0.57) fm for the copper (carbon) target. 

We modified the resonance shape of the (f> meson with 
the parameters ki and , then fitted again the observed 
mass spectra with the same procedure as before. For 
/cf^ we examined the two cases as described above. The 
best-fit results in both cases are shown in Fig. [31 The 
modification parameters ki, /c2°* and are common to 
the six spectra. Figure S] shows the confidence ellipsoids 
for the variation of with fci and fcj"* in the case (i) 
and (ii). In the case (i), obtained best-fit parameters are 
ki ^ 0.034j;[!™?, kl°^ = 2.6ti;^, and the minimum 
(Xmin) is 316.4. In the case (ii), best-fit parameters are 
ki = OmStoml, k\°^ = 0+5■^ and the x^i„ is 320.8. In 
both cases the Xmin ^as obtained with parameter ki ~ 
0.034, meaning the pole mass of the (f) meson decreases 
by 3.4 % at normal nuclear density. The x^j^ in the case 
(ii) (=320.8) is larger than that in the case (i) (=316.4) 
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FIG. 3: The best fit results in the case (i) (solid line) and 
(ii) (dotted line). The shift parameter fci, and broadening 
parameters ^2°' and A;|*, are common to the six spectra. The 
values of x^s are the sum of the squares of the deviations over 
54 data points in each panel. 
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FIG. 4; Confidence ellipsoids for the modification parameters 
ki and ^2°' in the case (i) in (a) and (ii) in (b). The values 
of Ax^'s in both panels are the differences from the Xmin(— 
316.4) at the best fit point in the case (i) which is shown 
by the cross in the panel (a). The best fit point in case (ii) 
is shown by the closed circle in the panel (b), and also in 
(a) since the ordinates are common to the both cases in the 
parameter space. 



by 4.4. When we fitted only the Cu data in the lowest 
(3"f region, where the major discrepancy occurs, the case 
(ii) was rejected at 99% C.L. and the best-fit parameters 
of fci and ^2°' do not change, while the case (i) was not 



In summary, we investigated the mass modification of 
the (j) mesons by studying the e^e^ invariant mass dis- 
tributions obtained in 12 GeV p + A reactions. The data 
obtained with a copper target revealed a significant ex- 
cess on the low-mass side of the (j> nieson peak in the 
(3j^ < 1.25 region. This observation is consistent with 
the picture of the (p modification in the nucleus. Thus, 
we conclude that in addition to our earlier publication on 
p/u modification this study experimentally verified 

vector meson mass modification at normal nuclear den- 
sity. 
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rejected. The data thus favor an increase of 
factor of 3.6 at normal nuclear density. 
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